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Abstract 

A fundamental task in developmental biology is to identify the mechanisms which drive mor- 
phogenesis. In many cases, pattern formation is driven by the positional information determined 
by both the gradient of maternal factors and hard-wired mechanisms embedded in the genome. 
Alternative mechanisms of positional information that contribute to patterning are the influence of 
signals derived from surrounding tissues. In this paper, we show that the interplay of geometrical 
constraints imposed by tissue shapes and hard-wired mechanisms, computationally implemented 
by a gene network model, can able to induce stable complex patterns. The rise of these patterns 
depends strongly on the geometrical constraints such as the shape of the surrounding tissues. 
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I. INTRODUCTION 



Developmental biology has identified several organizing principles that contribute in an 
orchestrated manner to embryogenesis and organogenesis. Several cell communication pro- 
cesses, such as lateral inhibition embryonic induction Q], cell growth and death ^| and 

organs. Such the processes are con- 
trolled by the information contained in the genome, which is usually the same in all cells. 
Therefore, a question arises: how do complex patterns emerge from initially homogeneous 
conditions? One way to address this issue is by morphogenic gradients. Morphogens are 
signaling molecules that can induce several distinct cell fates in a concentration-dependent 
manner P, |(|. The diffusion of a morphogen from a localized source generates a graded 
distribution across a field of cells. The local concentration of this molecule would induce the 
fate of nearby cells in a position-dependent manner. Typically, a morphogen is produced 
by a distinct localized set of cells, from which it moves into surrounding tissues. If cer- 
tain conditions are met, a spatial gradient of morphogen concentration, decaying away from 
the source, will result. When a signal is bound to its receptor, a specific intracellular sig- 
nal transduction pathway is triggered, leading to transcriptional and/or pos-transcriptional 
changes in responsive cells. If cells can sense the local concentration accurately, then spa- 
tial patterns of cell response can be generated. For example, in Xenopus embryos, it has 
been established that bone morphogenetic protein controls neural and epidermal fates in the 
vertebrate ectoderm, under the control of antagonists secreted by the organizer region of 
the mesoderm Q. Also wingless, hedgehog, and decapentaplegic proteins in the Drosophila 
wing imaginal discs, are well studied examples of secreted signaling proteins 8|]. 

From a mathematical point of view, Alan Turing showed a pathway to pattern formation 
without pre-patterns for its initiation P, Q]- This self-organized mechanism can be thought 
of as a competition between activation by a slow diffusing chemical and inhibition by a faster 
chemical, and has been largely applied to explain biological pattern formation (for example 
see P,lU(. However, it is important to keep in mind that biological strategies of patterning 
can be very different from the self-organized Turing-type mechanisms. Segmentation in 



Drosophila 3 and somitogenesis 14, [3| 



in vertebrate seem to follow a hard-wired 



strategy hierarchically organized from small regulatory gene networks. The extent to which 
patterns under development are self-organized or hard-wired remains an issue of debate 
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lalialial. 

In this paper we will explore, by means of computational simulations, the role of chemical 
signaling derived from surrounding tissues to patterning processes. In particular, we will 
show that specific geometrical constraints, like the shape of morphogen sources, are able 
to generate nontrivial endogenous gradients which will play a key role in the formation of 
complex structures. Modeling pattern formation is a current research topic in developmental 



biology 



20, 



2l| and to the best of our knowledge, the interplay of geometrical 



constraints, imposed by the shapes of surrounding tissues, and feedback control implemented 
in a gene network model has never been addressed before. Some studies about the effects 
of curvature and other geometrical constraint on pattern formation have been described in 



22, 



23]. 



In most cases, the shape of endogenous morphogen gradients is rarely known because the 
absolute concentration of the factor under analysis is so low. Furthermore, the regulatory 
gene networks, as well as the signalling pathway involved in developmental programs are 
not, in most cases, quantitatively characterized. For these reasons, in order to illustrate 
and to test the feasibility of our approach, we use a mathematical model and computational 
analysis. Our mathematical model sticks to the hard-wired strategy, implementing several 
negative feedback loops in a generic regulatory gene network. This gene network embedded 
in endogenous gradients is able to induce stable complex patterns. The rise of these patterns 
depends strongly on the geometrical constraints such as the shape of the surrounding tissues. 

Successful development can take place in a range of environments. Developmental mech- 
anism of patterning must be robust and precise. Signals produced at wrong place or time 
can lead to an inappropriate developmental process. The relevance of noise and robustness 
will be also discussed. 



II. METHODS 



A. Shaping morphogen gradients 



Although pattern formation is indeed induced by graded activation of a signaling path- 
way, several factors are involved in the distribution of a morphogen in a developmental field 
The mechanisms involved in gradient shaping is one of the main issues debated in 



way, 
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developmental biology today. Such mechanisms include passive diffusion in the extracellular 
matrix j^J, repeated rounds of endocytosis and exocytosis 26], cytonemes j^, argospmes 
28 1, and a recently proposed mechanism by mRNA decay in growing structures [29]. It 
is beyond the scope of this paper modeling the shape gradient genesis. We will adopt the 
simple localized source-dispersed sink models, where the shape of a morphogen gradient is 
dictated by the rate of diffusion away from the site of synthesis, along with the rate degra- 
dation. The gradient can be broadened by either increasing the diffusion rate or decreasing 
the degradation rate, while the opposite effects can produce a much steeper gradient. In 
particular, the diffusive substances are secreted by the fixed boundary of surrounding tis- 
sues. Thus, we also consider that these substances are degraded over all developmental field. 
We can represent these processes by the following equations 

a = D a V 2 a — X a a 

h = D h V 2 h-X h h, (1) 

where a (h) denotes the activator (inhibitor) concentration which is released by the tissue 

1 (tissue II) and degraded over the domain at rate A a (A/J. Fig. 1 shows a schematic 
representation of surrounding tissues, boundaries (bottom panel) and the concentration 
profiles a and h over a horizontal segment in the middle of a semi-elliptic domain (top 
panel) . 

In order to evidence the key role of geometrical constraints in the patterning process, the 
overall processes were implemented on both semi-elliptical and semicircular domains. Fig. 

2 depicts the concentration fields of activators (A) and inhibitors (B) in a portion of the 
semi-elliptical domain. For the semicircular domain the concentration fields of activators 
and inhibitors are shown in (A) and (B) respectively and they were obtained using the 
parameters in Table 1. Despite of small differences between the top and bottom morphogen 
distributions, they are sufficient to induce completely different patterns of expression of 
regulated genes, as we will show later. 

B. Hard-wired model for development 

The striped pattern of pair-rule genes in Drosophila is a consequence of the fact that 
each stripe is separately controlled by dedicated transcription regulation [l^. This is an 
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example where biological strategies of patterning seem to follow a hierarchically organized 
building strategy from small modular regulatory gene networks. We will construct a small 
regulatory network (see Fig. 3) which is able to produce complex patterns from an ini- 
tially homogeneous distribution of the genes products. Our mathematical model considers 
several biochemical and biophysical processes currently present in many biological descrip- 
tions, including: diffusion, degradation, noise and of course, regulatory activation and/or 



feedback loops, which seem 
. We will describe reg- 



repression of gene expressions. The last one forms two negative 
to be essential control mechanisms in developmental networks 
ulatory activation sigmoidal functions s + (x,8,n) = x n /(x n + 9 n ). If the concentration of 
activator factor is below the threshold 9, the gene is poorly expressed, whereas above this 
threshold its expression grow until saturate. The regulatory inhibition will be represented 



by s = 6 n /(x n + 6 n ) = 1 — s + . Sigmoidal functions, in 



20 



21 



particular Hill functions, have 



>'l, Eill • I ' general, there are 



been used in a variety of genetic regulatory models 
several regulatory interactions acting over a single gene, and in these cases the resulting 
expression rates will be described by the product of sigmoidal functions corresponding to 
each interaction [^, rather than to the additive formulation used in 21, For the 
sake of simplicity, our model does not consider other factors as growth, cell division, pos- 
transcriptional regulation, etc., and does not explicitly distinguish between genes and their 
products. 

We will consider five genes which form a small network of interacting genes. Their product 
concentrations, represented by U{ with i = 1, ... ,5, evolve following the reaction-diffusion 
equations 

Ui = diV 2 Ui + fi (a, h, {ui}) - X Ui Ui + rje (t) u h (2) 

where fi represents the regulatory transcription control of gene i, and Aj is the degradation 
rate. The last term represents a multiplicative noise, e is a random variable that is assumed 
uncorrelated and Gaussian distributed (null mean and variance unitary), t] represents the 
intensity of noise. Eq. (2) governs the spatio-temporal evolution of the small regulatory 
gene network represented in Fig. 3. This network suffers the influence of the surrounding 
tissues through the biochemical agents concentration of a and h (which are constant at a 
time). The particular form of fi depends on the regulatory inputs acting over gene i, and 
will be discussed later. 

In the hierarchically organized hard-wired strategy, the regulatory role of the genes be- 



longing to a network could be organized in different levels, for example: positional level, 
regulatory level, output level, etc.. These levels interact between them to constitute a spe- 
cific network structure. In this framework, we can understand a complete developmental 
program as a network of interacting gene networks. 

In this spirit, the positional level will be assigned to Ul, the gene associated to concentra- 
tion u\. Ul will decode the positional information derived from the geometrical constraints 
(boundaries) through a and h concentrations. The transcription of gene Ul will be activated 
when a reaches a threshold 6> 1(l and will be inhibited by h and its own product at the thresh- 
olds 9 lh and #n respectively (the top panel of Fig. 1 shows the values of these parameters). 
Therefore, the expression of ui is governed by Eq. 3: 

/i(a, h,ui} = n s + (a,9 la ,n la ) s~ (h,9 lh ,n lh ) x 

s~(ui,0ii,nn). (3) 

The expression of Ul activates in a concentration-dependent manner the genes U2, U3, 
and U4, (associated to the concentrations u 2 , u 3 and it 4 respectively), which constitutes 
essentially a regulatory level of the hierarchy. U2, U3, and U4 are activated at thresholds 
#21 > #31 > #41, respectively, and each product inhibits the transcription of the other genes, 
thus forming two negative feedback loops. So, we can write: 

h ({ui}) = r 2 2 s + (u 2 , #22, n 22 ) + r 2 i s + (ui, 9 2 i,n 2 i) x 

S~ (U 3 , #23, n 23 ) S~ (u 4 , #24, n 24 ) , (4) 

where {u^} represents the set of variables {u\, u 2 , u 3 , 114, u 5 }. for r 2 2 > this gene is also 
self activated. We have added a self activated term in (4) to study their effects over noise 
control. 

/3 ({«<}) = r 3 S + (ill, 3 1, r^l) S~ (M 2 ,#32,n 32 ) X 

S~ (M 4 ,#34,n 34 ) , (5) 

U d U i}) = r 4 S + (til, #41, n 4 l) S_ ( M 2, #42, n 42 ) X 

(«3, #43,^43) • (6) 

Thus the gradient associated to u x generates discrete spatial domains of regulatory factors 
which interact between them and generate a more complicated spatial domain of transcrip- 
tion of the third level of the hierarchy, the output level. For simplicity, we consider also 
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that this level is constituted by only one gene. Its activation is promoted by u\ and it is 
down-regulated by the gene products of the former level, u 2 , M3 and u±. 

h {{ u i}) = r 5 s + (ui, 6»5i, n 51 ) s~ (u 2 , 052, n 52 ) x 



The positive and negative feedback loops of the model are shown in Fig. 3. It is important 
to keep in mind that the model was constructed to illustrate that gene networks, embedded 
in feedback control principles, can induce or not complex patterns depending on geometrical 
constraints of the problem. For the aim of our study, this small network is enough. Of 
course, other kinds of network can be designed to model specific developmental systems. 

III. RESULTS 

A. Geometrical constraints and gene regulatory network 

The boundary problem is completely specified by Eqs. (1), the initial condition (IC) and 
the boundary condition (BC). In this sense, a = 1 over C and null over C", the two vertical 
boundary segments are connected by a periodical boundary condition, while the inhibitor 
is h — 1 over C and null over C, the two vertical boundary segments are connected by 
periodic boundary condition. Initially, both a and h concentrations are null over all domain 
compatible with BC. We integrate numerically the Eqs. (1) for two different situations 
which differ in degradation parameter values. The simulations run for a long period of 
time (5,000 steps) after which the resulting concentration a and h are stationary. Then, 
these concentrations will be used as input signals acting over the network in several different 
conditions. 

The aim of the paper is to study the effects of geometrical constraints on the patterning 
in developmental motivated models rather than to modeling a particular system. For this 
reason it is not necessary to fit the parameters to reproduce experimental observation. How- 
ever, as the particular model chosen for our purpose has 52 parameters, we need to establish 
some guidelines so as to reach stable stationary patterns before implementing massive com- 
putation. In this model, the gene expression rates involve Hill functions whose coefficients 
riij range from 1 to 5. n^j > 1 reflects the fact that ds-regulatory systems have usually 




(7) 
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many binding sites for each trans-regulatory element. These parameters can influence the 
shape patterns, in particular by controling the sharpness of the patterns. The 0ij parameters 
represent the threshold for switching on or off the transcription processes. They were setted 
in the interval [0.0,1.2]. Some patterns are very sensitive to these parameters while others 
are robust. In particular, the thresholds corresponding to the activation of U2, U3, and U4 
by Ul are sorted in decreasing order so as to form a sequence of activation; we setted them to 
#4i = 0.4, #3i = 0.6 and 6> 2 i ranges from 0.9 to 1.2. The diffusion constants and degradation 
rates are initially the same for all biochemical species (0.01 and 0.20 respectively). Initially, 
the level of noise was null for all species. Table 1 and 2, display the parameters values used 
for most computations, otherwise we will mention other parameters values. 

The set of partial differential equations were discretized using standard finite difference 
methods. The resulting large-scale dynamical system was integrated in time using forward 
Euler integration. We have performed simulations for different sets of parameters subject 
to the above mentioned restrictions. We have found that patterns are sensitive to the 
spatial profile concentration of a and h and the parameters that regulate the activation 
and repression of Ul expression, 9 ia , Q^h and 0u. For any given profiles a and h, there 
are regions in the parameters space where no patterns form, mainly due to the fact that 
activator concentration a does not reach the threshold 6\ a . An interplay of these parameters 
and both a and h profiles control the position and the coarse-grain aspect of the pattern 
shapes. 

Figure 4 depicts the density of concentrations of u 2 (A), u 3 (B), w 4 (C) and u 5 (D) 
obtained using the input profile 1 (see Table 1) after 400 time steps, when stationarity is 
reached. We can see that the concentrations present complementary aspects. In general, the 
pattern of the output level is particularly complex as it is under the effects of four regulatory 
fields. In contrast, the concentration corresponding to u 2 is the simplest one because the 
regulatory field is almost homogenous in its expression domain. We have observed trough 
several simulations that these characteristics are quite general. 

The remarkable observation in our study is the fact that geometrical constraints play a key 
role when surrounding tissues secret morphogenic substances. Localized patterns, similar to 
those developed in the semi-elliptic domain, were not formed in a semi-circular domain. In 
the last case, we were not able to tune the parameters in order to obtain patterns different 
from semi-circular stripes as shown in Figure 5. The reason for that is based on the fact that 
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(despite of diffusively substance a and h are uniformly secreted by the surrounding tissue) 
the curvature of the tissues is not homogenous in the elliptic case, but it is homogenous 
in the circular case. The two fronts of secreted substances are able to generate a localized 
expression domain for Ul in the non-homogenous curvature case. The position where the 
complex structures develop in Figure 4 corresponds to the maximal value of the surrounding 
tissue curvature. In contrast, for homogenous curvature the expression domains are not 
localized and that does not depend on tuning adequately the network parameters. 

The same network can generate different patterns depending on the parameters values 
and the input profiles. Fig. 6 illustrates four patterns of expression level corresponding to 
W5 after 400 steps. The patterns displayed at the top (A and B panels) have been obtained 
using profiles 1 (see Table 1) and two different sets of parameter values. In panel A the 
system reached stationarity, while panel B reached stationarity after 800 steps (see Movie 
1). Stationary patterns displayed at the bottom (C and D panels) have been obtained using 
profiles 2 (see Table 1). The parameter values are shown in Table 3. Often different sets of 
parameter values lead to the same expression pattern, but the way to reach this pattern can 
differ from set to set. The movie 1 illustrates the evolution of the network in two different 
conditions. The right panel of the movie corresponds to the same condition as in Fig. 6B, 
while in the left panel we changed slightly the value of parameter to r% = 32.1 x 10~ 2 . Both 
simulations correspond to 800 steps. We can see that the final patterns are almost the same, 
but their evolutions are quite different. Small perturbations of boundary conditions affect 
the patterning processes differently, depending on the network parameters. For example, 
Fig. 7 shows patterns generated by the network in the same conditions as in Fig. 6, but 
with a perturbation in the boundary tissues (both in C and C). The panels of Fig. 7 
illustrate that the perturbation can affect a little, as in panel A, or too much, destroying 
completely the patterns as in panels C and D. This result suggests that the shape of the tissue 
strongly influences the formation of pattern expression and biochemical interactions are not 
a necessary condition to explain alterations in the resulting gene expression patterns. The 
influence could be through the geometrical constraint imposed by the surrounding tissues. 

In general, genetic network which underlie the developmental programs of living cells, 
must withstand considerable random perturbations. This occurs as fluctuation in, for ex- 
ample, transcription, translation, and RNA and protein degradation 



33, 



We have also 



studied by numerical simulation the effects of the noise on the patterns in two situations. For 
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all the previous cases we have setted r 22 = 0.0, i.e. there is no autocatalytic feedback acting 
over w 2 . It would be interesting to compare the effect of noise when this positive feedback 
loop is present in the network. Figures 8 and 9 show the expression profiles of u% and u$ re- 
spectively, obtained for five different levels of noise rj = 0.00 (A and F), rj — 0.05 (B and G), 
7] = 0.10 (C and H), 77 = 0.15 (D and I) and rj = 20 (E and J). In both figures the top panels 
(A-E) correspond to the network without the above mentioned autocatalytic loop, while the 
bottom panels correspond to the network with the loop. The network parameter values used 
in both cases were the same as those used in Figure 7A, except for parameters r 21 and r 22 , 
which in the case of the network with the loop, change from 6.85 x 10~ 2 to 4.84 x 10~ 2 , 
and from zero to 4 x 10~ 2 respectively. We can see that pattern structure remains almost 
unaltered for noise level smaller than ~ 10% for the network with the autocatalytic loop. 
On the other hand, when this loop is absent (r 22 = 0.0), the patterns seem to be more 
noise sensitive. This example illustrate that the network is robust against a considerably 
high level of noise. This robustness could arise from the fact that there are two negative 
loops underlying the regulatory interactions. Furthermore, this example shows, as in many 
homeotic genes, that the positive feedback loop helps the maintenance of their expression, 
thus contributing to the stabilization of the patterns. 



B. Geometrical constraints and self-organized strategy 

In order to study the influence of geometrical constraints on a self-organized system type, 
we have also implemented the Gray-Scott model on the semi-circular and semi-elliptical 
domains. The Gray-Scott model is a variant of the autocatalytic Selkov model of glycolysis 
which corresponds to the following two irreversible reactions: 

U + 2V -> W 
V -> P. 

P is an inert product. Both U and V are removed by the feeding process. A variety of 
spatio-temporal patterns were derived in response to finite-amplitude perturbation [35]. We 
consider that u and v are coupled to the activator field a and an to inhibitor h as follows: 

u = D u V 2 u - uv 2 + F(l -u)+r u f (a, h) 

v = D v V 2 v + uv 2 - (F + k)v + r v g (a, h) , (8) 
10 



where 



f(a,h) = s + (a, 9 u>a , nu t a) s (h,9 Uth ,n Uth ) 
g(a,h) = s + (a,9 V:a ,n V:a ) s~ (h,9 Vjh ,n v>h ) . 

We study numerically the response to different positioned perturbations of the system 
which evolve in both geometries mentioned above. The initial condition corresponds to the 
trivial state (u — 1 and v — 0). A 16 by 16 mesh point area located over the semimajor 
(horizontal) axis was perturbed to u — 0.50 and v = 0.25. The system suffers the influence 
of a and h fields determined by the profile 1 (see Table 1). Figure 10 shows the patterns 
obtained from the system described by Eqs. (8) on the elliptic geometry (for system response 
on the circular geometry see Fig. 11). The horizontal position of perturbations were centered 
at 38.3, 42.3, 46.3, 50.3 and 54.3; while parameter k ranges from 0.060 to 0.066 and is 
uniformly spaced. Other adopted parameters are shown in Table 4. As we can see, for the 
several values of the parameter k, the final patterns upon 2,000 time steps depend strongly 
on the position where the perturbation occurs; while the geometry of the domain, either 
elliptical or circular, seems to have very little influence. It should be mentioned that, in 
contrast to the patterns derived in the previous section, these patterns are not stable, they 
continue growing until all the domain is filled. 



IV. DISCUSSION 



Early development of multicellular organisms is marked by a rapid initial increase in their 
cell numbers, accompanied by morphogenetic processes leading to the gradual formation of 
organs of characteristic shapes. During morphogenesis, through differentiation under strict 
genetic control, cells become more and more specialized. Further, genetic mechanisms as 
morphogenesis also require generic physical principles, such as, diffusion, spreading, dif- 
ferential adhesion, chemotaxis, etc.. As a consequence, development rely on an intricate 
interplay of generic and genetic mechanisms. In this paper, we have addressed the issue 
on how chemical signaling derived from surrounding tissues can drive patterning processes. 
Our results, derived by computer simulation, suggest that the shape of the source (tissue) 
and other geometrical constraints strongly influence the formation of complex structures. 
In particular, we found that non-homogenous curvature of signaling tissues can generate 
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complex patterns when acting on gene network, while the same network embedded in a 
geometry with homogenous curvature (boxes and circles) does not form localized structures. 
This result suggests an important consequence for development: let us consider that the 
shape of an organizer (tissue I or II in our case) is controlled by a gene network A, and the 
genie response to the organization field is controlled by a network B (the network shown 
in Fig. 3 in our case), which is not regulated by genes of network A. Any mutation in A 
that alters the shape of the organizer and consequently the organization field, will affect the 
development of the induced structure or organ. Thus, we can conclude that biochemical 
interaction is not a necessary condition to explain alteration in the overall induced output; 
the interaction could be underlaid by a geometrical constraint imposed by the organizer. 
In contrast to this hard-wire system, the self-organized system examined here (Gray-Scott 
model) seems to be more insensitive to the gradients derived from surrounding tissues. 

We have also studied the robustness of patterns in response to noise. In this sense, two 
kinds of network have been used to illustrate the effect of noise in either the presence or 
absence of a positive feedback loop. The pattern driven by the first network seems to present 
higher robustness to noise than the second one. However, the two negative feedback loops 
presented in both cases were able to guarantee the stability of the pattern structure against 
a considerably high level of noise. 

Although we have used a particular network topology with several mutually-inhibiting 
factors, there are many other possibilities for the topology that could produce interesting 
patterns. Each topology can response in a different way to the intrinsic noise of the pattern- 
ing process. The particular case examined here was enough to conclude that the hard wired 
mechanism, as the regulatory networks can create complex biological shapes out of a simple 
structure. This mechanism could be enriched by incorporating geometrical constraint as key 
ingredients of morphogenesis processes. 
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Tables 



TABLE I: Parameter values for diffusion, degradation and geometrical constraints used to obtain 
the input profiles 1 and 2. These profiles differ in the degradation values 7 a and jh which were 
reduced to 0.5 x 1CT 3 to generate the input profile 2. 



Parameters 


Values 


Parameters 


Values 


di 


0.01 


Semi-ellipse 




H (* + 5) 


0.20 


internal semimajor axis 


40.00 


75 


0.35 


internal semiminor axis 


14.14 


D a 


0.04 


external semimajor axis 


86.60 


D h 


0.04 


external semiminor axis 


31.62 


la 


1.0 x 10~ 3 


Semi-circular 




lh 


1.0 x 10~ 3 


internal radius 


24.70 


At 


1.00 


external radius 


51.96 


Ax 


0.40 







15 



TABLE II: Parameter values of the gene network used to obtain the pattern shown in Figure 4. 
These values were also used in other simulations whenever we have not stated otherwise (see also 
Table III). 



Parameters 


Values 


Parameters 


Values 


n 


6.25 


r 4 


5.8 x 10~ 2 


Ola 


0.40 


#41 


0.40 


Glh 


0.05 


#42 


0.25 


#11 


0.20 


#43 


0.20 


Ula 


1.00 


"41 


4.00 


nih 


3.00 


"42 


2.00 


nn 


5.00 


"43 


2.00 


T21 


8.06 x 10~ 2 




20/27 


021 


1.05 


#51 


0.55 


#23 


0.35 


#52 


0.20 


#24 


0.45 


#53 


0.30 


nil 


4.00 


#54 


0.50 


"23 


2.00 


"51 


4.00 


"24 


2.00 


"52 


3.00 


T3 


0.21 


"53 


3.00 


031 


0.60 


"54 


3.00 


032 


0.10 


PFL parameters 


#34 


0.10 


^22 


0.00 


"31 


4.00 


#22 


0.10 


"32 


2.00 


"22 


3.00 


"34 


2.00 
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TABLE III: Parameter values of the gene network used to obtain the patterns shown in Figure 7. 
The values of parameters not shown in this table are the same as in Table II. 



Fig 


. 7 A 


Fif 


5. 7B 


Fi£ 


r. 7C 


Fig 


. 7D 




Parameters 


Values 


Parameters 


Values 


Parameters 


Values 


Parameters 


Values 




6.85 x 1(T 2 


n 


9.25 


n 


7.50 


n 


7.50 x 


10~ 6 


021 


0.99 


0n 


0.265 


01a 


0.20 




8.46 x 


10-2 


"21 


5.00 


T21 


12.09 x 10~ 2 


T21 


12.09 x 10^2 


021 


0.95 






0.215 


021 


0.85 


021 


1.30 


"21 


5.00 




031 


0.65 


"21 


5.00 


"21 


5.00 


^3 


0.27 




r 4 


5.60 x 10~ 2 


^3 


32.02 x lO- 2 


^3 


0.18 


"31 


2.00 






0.18 


r 4 


6.24 x 10-2 


r 4 


5.20 x 10-2 


r 4 


9.75 x 


10-2 


051 


0.50 




35/216 




10/9 


041 


0.45 




052 


0.55 


051 


0.50 


051 


0.52 


042 


0.20 




053 


0.50 


052 


0.40 


052 


0.40 


"41 


2.00 




054 


0.20 


"51 


2.00 


053 


0.30 


T5 


8.54 




n 5 i 


2.00 






054 


0.25 


051 


0.50 












"51 


3.00 


052 


0.40 
















053 


0.20 
















054 


0.20 
















"51 


2.00 
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TABLE IV: Parameter values for the Gray-Scott model. 



Parameters 


Values 


Parameters 


Values 


F 


0.04 


C(= Ouh) 


0.50 


r u 


0.001 


@va = @vh) 


0.50 


r v 


0.002 


n ua{ = n uh) 


1.0 


D u 


0.0336 


n va{ = n vh) 


1.0 


D v 


0.0168 
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Figures and movie 




FIG. 1: A schematic representation of surrounding tissues, boundaries (bottom) and the concentra- 
tion profile of activator a and inhibitor h over a horizontal segment in the middle of a semi-elliptic 
domain (top). 
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FIG. 2: 3D plots of the concentration of a, secreted by the tissue I (A and C), and h (B and D) 
which is secreted by the tissue II, obtained from Eqs. 1, after 5,000 time steps. The panels A and 
B correspond to the elliptical geometry, and panels C and D correspond to the circular geometry. 
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FIG. 3: Schematic diagram of the network. Circles indicate genes, blunt-end arrows denote inhi- 
bition and pointed arrows denote activation. Regulatory activation is represented by arrows and 
regulatory inhibition by small perpendicular edges, as usually. This network presents two negative 
feedback loops, one where U2 inhibits U3 that inhibits U4 that inhibits U2; and another which is 
formed by U4 that inhibits U3 that inhibits U2 that inhibits U4. There is also an auto-inhibitory 
loop acting over Ul and an autocatalytic loop acting over U2. 
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FIG. 4: Density plot of concentrations of U2 (A), u% (B), (C) and u§ (D) after 400 steps in 
a semi-elliptic domain. The model parameters used to compute these patterns are displayed in 
Table. 
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FIG. 5: Density plot of concentrations of U2 (A), U3 (B), U4 (C) and U5 (D) after 400 time steps. 
In this semi-circular domain, no localized structure were observed for any set of parameter values. 
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FIG. 6: Density plot of concentrations of u§ after 400 obtained for inputs profile 1 (A and B) and 
inputs profile 2(C and D) with the network parameters displayed in Table 3. 



24 



A _J 


? ] 


3|] 

c _J 


D _J 



FIG. 7: Density plot of concentrations of u§ after 400 steps obtained in the same condition that 
Fig. 6 with the exception in the boundary condition which were perturbed. 
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FIG. 8: Density plot of 113 obtained, after 400 steps, when the dynamics of the gene network is 
perturbed with noise. The top panels correspond to patterns obtained from the network without 
the positive feedback loop for increasing level of noise (sorted from left to right). The bottom panels 
correspond to patterns obtained from the network with the feedback loop for increasing level of 
noise (sorted from left to right). The level of noise used were 77 = 0.0 (A and F) 77 = 0.05 (B and 
G), 77 = 0.10 (C and H), 77 = 0.15 (D and I), and 77 = 0.20 (E and J). For further information, see 
the text. 
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FIG. 9: Density plot corresponding to U5 obtained in the same situation as in Fig. 8. 
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FIG. 10: Elliptic Geometry: Density plot of the v concentration for different position of perturba- 
tion, increasing from left to right, and for different values of parameter k which increases from top 
to bottom (see details in text). 
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FIG. 11: Circular Geometry: Density plot of the v concentration for different positions of pertur- 
bation (increasing from left to right) and for different values of parameter k (increasing from top 
to bottom). 
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Movie 1: (Please see at http:/ /glia.if.sc.usp.br/luis/teste.mov| The left panel corresponds to 
the temporal evolution of the pattern Fig. 7B during 800 steps. The right panel corresponds 
to the evolution in the same condition with the exception of pattern r% which changes from 
32.02 x 10~ 2 to 32.1 x 10~ 2 
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